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ABSTRACT 

Photometric redshifts are used to determine the rest frame luminosity function (LF) 
of both early-type and late-type galaxies to Mb ~ -17.6 for the cluster C11601-I-42 
at z = 0.54. The total LF shows a steep faint-end slope a ~ -1-4, indicating the 
existence of a numerous population of dwarf galaxies. Luminous galaxies, with Mb ^ — 
19.5 are mostly red, early-type galaxies, with a LF best described by a Gaussian. 
Faint galaxies are predominantly blue, late-type galaxies, well fitted by a Schechter 
function with a ^ -1.7. Compared to clusters at lower redshift, the steepening of 
the faint end starts at brighter magnitudes for C11601-I-42, which may indicate a 
brightening of todays dwarf population relative to the giant population with increasing 
redshift. Early-type galaxies are centrally concentrated, and dominate the core region, 
implying that the radial gradient of early-type galaxies seen in local clusters is already 
established at z 0.5. Bright, late-type galaxies are rare, consistent with a decrease 
in star formation in field galaxies as they are accreted on to the cluster, while faint, 
blue galaxies are evenly distributed across the cluster, except for a depletion in the core 
region. The blue fraction is /s ~ 0.15, which is somewhat lower than the Butcher- 
Oemler average at z ^ 0.5. The value of Jb is found to increase with limiting 
magnitude and with radius from the centre. 

Key words: Galaxies: clusters: individual: C11601-I-42 - galaxies: distances and red- 
shift - galaxies: luminosity function, mass function - Cosmology: observations 



1 INTRODUCTION 

Butcher & Oemler (1978a; 1984, hereafter B084) found that 
the fraction of blue galaxies increases dramatically between 
redshifts 2 = and z ~ 0.5, now usually referred to as the 
Butcher-Oemler effect (hereafter BO effect). This t rend has 
been confi rmed by e.g. Rakos & Sch ombe rt (1995), EUing- 
son et al. ( 2001), Kodama & Bower (2001) and Margoniner 
et al. (2001). The increase of blue galaxies is explained by 
a higher rate of star formation in cluster galaxies at earlier 
epochs. The cause of this higher star formation rate is, how- 
ever, not fully understood. Butcher & Oemler (1978a) noted 
that the excess of blue galaxies at higher redshifts compared 
to todays clusters approximately matches the excess of SO 
galaxies in local clusters. If such a change in the morpholog- 
ical mix is caused by an evolution of the blue galaxies into 
SO galaxies, once the gas in the star forming blue galaxies is 
exhausted or tidally stripped, this could lead to the observed 
BO effect. The evolution of the fraction of SO's between high 
and lo w redshift clusters is dis cussed in e.g. Dressier et al. 
( [L997| ) and Fasan o et al. ( ^000| ). 



Bower (1991) and Kauffmann (1995) explain the BO 
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effect as a natural consequence of the hierarchical cluster- 
ing model for structure evolution in the Universe. In this 
picture structures, such as clusters, are built up from bot- 
tom up, via repeated mergers of subunits. Bower shows that 
the infall rate of field galaxies to clusters increases with look- 
back time, leading to a higher fraction of blue galaxies in the 
clusters. The trend found is not as strong as in the B084 
observations, but can be fitted if the blue fraction in the 
field increases with redshift. 

Kauffmann shows that high redshift clusters assemble 
during a much shorter time interval than comparable clus- 
ters today. This in turn implies a higher rate of infall of 
galaxies from the field at higher redshift, producing the 
BO trend. However, the strong increase in the blue fraction 
of the observed clusters is explained only in high-density 
(r^M = 1) CDM models, or MDM models. Low-density 
CDM models have difficulty in explaining the strong trend 
at 2: < 0.4 ( Kauffmann 1995 ). Recent simulations by Di- 
aferio et al. (2001) find that a BO trend is found also in a 
low-density model (57 Af = 0.3, JIa = 0.7), but significantly 
weaker than in observations. 

It is not clear if star formation is triggered in the in- 
falling galaxies, possibly creating a burst, or if star forma- 
tion in the field galaxies is simply quenched over some period 
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of time. Both these scenarios seem compatible with the ob- 
served BO effect. 

The resuhs from B084 also show that open clusters 
have a larger fraction of blue galaxies than the more cen- 
trally concentrated, compact clusters. (A formal definition 
of open and compact clusters is given in Section 4.1.) This 
can be explained in the context of the hierarchical scenario, 
if the smaller concentration in open clusters is a reflection 
of a more recent assembly of its members. 

The presence of substructure in a cluster can be a sign of 
recent subcluster-subcluster merger. From a study of three 
clusters with signs of s ubstr ucture at 0.1 < z < 0.2, Mete- 
vier, R omer fc Ulmer (|200C|) find that at least two of these 



have ai exceptionally large fraction of blue galaxies. This 



supports a hierarchical scenario where mergers induce star 
formation in the cluster galaxies. 



Furthermore, results from Margoniner et al. ( 2001) indi- 
cate that poor clusters in general have higher blue fractions 
than rich clusters at comparable redshifts. It is therefore 
clear that the BO effect is not a simple redshift effect, but 
must be viewed as a consequence of several factors. These 
include richness, central concentration, epoch of formation 
and presence of substructure. 

A further prediction of the hierarchical CDM model is 
that the number density of galaxies should be dominated by 
low mass systems. This is reflected in the luminosity function 
(hereafter LF), defined as the number density of galaxies per 
unit luminosity, or per unit magnitude. For the hierarchical 
scenario, the LF is predicted to have a steep faint-end slope, 
refiecting the incr easing number of low mass systems with 
faint magnitudes (White & Frenk 1991). Therefore, studies 
of the LF, both in the field and in clusters, is one of the 
most direct probes for testing the hierarchical model. 

Detailed studies of the spatial and morphological distri- 
bution of cluster galaxies are important for understanding 
the internal properties of clusters, as well as the evolution of 
these over time. To do this observationally, it is necessary to 
isolate the galaxies in a given field belonging to the cluster, 
from the foreground or background field galaxies, projected 
on to the cluster. Spectroscopy is costly, since a cluster field 
may consist of thousands of galaxies. More important, for 
clusters at higher redshifts spectroscopy is limited to only 
the most luminous galaxies. The traditional way of correct- 
ing for field contamination at faint magnitudes is to use a 



cluster- Free field, preferentially close to the cluster, and cal- 
culate the number of galaxies as a function of magnitude 
to be subtracted from the cluster image. This method is, 
however, increasingly uncertain at higher redshifts, where 
the ratio of field galaxies to cluster galaxies increases (e.g. 
Driver at al. 1998a). In Section 3.7 we show that the field 
contamination can exceed 80 per cent of the total number of 
galaxies for clusters at redshifts z ^ 0.5. A further important 
disadvantage is that it is not possible to assign membership 
to a specific galaxy with any confidence. 

An alternative and more powerful way to determine 
cluster membership is to use photometric redshifts based on 
broad band photometry to estimate redshifts. The main ad- 
vantage with this method is that a large number of redshifts 
can be determined in a relatively short time, and to magni- 
tudes considerably fainter than traditional spectroscopy. An 
important bonus is that the photometric redshift technique 
also provides an estimate of the Hubble type for the ob- 



served galaxies, or more precisely, their broad spectroscopic 
characteristics. This information can be used, together with 
redshift information, to determine the fraction, as well as 
the spatial distribution, of different types of galaxies. We 
later discuss this technique in more detail. 

The trade-off with this method is that uncertainties 
in the estimated redshifts are considerably larger compared 
to spectroscopy. The accuracy has, however, improved and 
typical rms deviations between spectroscopic redshifts and 
broad band phot ometric redshifts are Azrma ^ 0.05 at 



0.6 < z < 1.0 (IConnoUy et al. 1995| ; |Liu fc Green 199^ 
Lubin fc Brunner 19991) and Azrms/(i + z) = 0.0 6 -0.1 up 



to z ~ 6 ( Fernandez-Soto, Lanzetta & Yahil 1999 



Benitez 



200C). This accuracy is sufficient for determination of clus- 
ter membership and absolute magnitudes, but is clearly not 
useful for dynamic studies. 

In this paper we demonstrate the use of photometric 
redshifts for a population study of the cluster C11601-I-42. 
With information on cluster membership and spectroscopic 
galaxy type, we calculate the luminosity function separately 
for early-type and late-type galaxies within the cluster, as 
well as their spatial distribution. Previously, photometric 
redshifts have been used to determine cluster membership 
for a few clusters only. The most extensive study of photo- 
metric redshifts appli ed to clusters of galaxies is presented 
by Brunner & Lubin (200C). They in paricular, discuss how 



the method can be used to either maximize completeness, 
or to min imize contamination from field galaxies. Lubin & 
Brunner ( [L999|) present results on the distribution of galaxy 
morphologies in two distant clusters, where cluster member- 
ship is det ermin ed by photometric redshifts. Kodama, Bell 
& Bower (1999) use photometric redshifts for determining 
cluster members in Abell 370 at z — 0.37. Connolly et 



al. (1996) include photometric redshifts when investigating 
a superclust er str ucture at z = 0.54, and Connolly, Szalay 
& Brunner ( 199^ ) use the method for measuring the evo- 
lution of the angular correlation function in the HDF. Fi- 
nally, Bezecourt et al. (2000) use photometric redshifts when 
analysing weak lensing of the cluster Abell 2219 at z — 0.22. 
For studies of the LF and population content of clusters this 
method has, however, only had limited use. 

The cluster C11601-I-42 is an intermediate rich cluster 
at z — 0.54. It was chosen because there already exis t both 
high resolu tion HST imaging of the core of the cluster ( 3mai] 



ot al. 1997), and a fairly large numbe r of spectroscopic red - 



shifts of galaxies in the cluster field (Dressier et al. 1999). 
The cluster is included in the MORPHS sample of ten clus- 
ters at 0.37 < 2 < 0.56 (SmaU et al. 1997). The available 
information is valuable both for testing the photometric red- 
shifts technique, as well as for determining the morphologi- 
cal properties of the cluster. 

The plan of this paper is as follows. In Section 2 we 
present the data, in Section 3 we describe the photometric 
redshift technique used. Results are presented in Section 4, 
followed by a discussion in Section 5. Finally, conclusions 
are given in Section 6. 

We assume in this paper a flat cosmology with Qm ~ 
0.3 and Q.a = 0.7, and a Hubble constant Hq — 50 km 
s~^ Mpc~^. The absolute magnitude for other choices of the 
Hubble constant are related by: M(j:/p=5o) = Mj^f^^gs) — 
0.57 = M(Ho=ioo) - 1-51. At the cluster redshift, z = 0.54, 
the absolute magnitudes for different cosmologies are re- 
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lated by: M(nj^j=o. 3,0^=0.7) = A^(nM=i.o,nA=o.o) - 0.42 



-^(nM=o.2,nA=o.o) —0.24. 



Magnitudes are given in the standard Vega based sys- 



tem. 



2 THE DATA 



2.1 The cluster C11601+42 

C11601+42 is located at R.A. 16''03'"09".8 and Dec. 
42° 45' 18"(J2000.0). The re dshift has been determine d 
spectroscopically to z = 0.54 (Oke, Gunn & Hoessel 1996), 
and the velocity dispersion is 1210 km s~^ dDressleT^t al 
199S). The cluster is a weak X-ray source with Lx (0.3-3.5 
0.35 /i"^ 



keV) 



10*" ergs (Henry et al. 1982; Small et al 



1997). The galactic extinction is cori sistent with Ab = 0.00 
(pchlegel, Finkbeiner fc Davis 1998|). 



2.2 Observations 

The cluster was observed with the 2.56m Nordic Op- 
tical Telescope (NOT) and the Andalucia Faint Object 
Spectrograph and Camera (ALFOSC) during four observa- 
tional runs 1997-1999 (Table ALFOSC contains a Lo- 
ral 2K^ chip with a field of 6'. 5x6'. 5 and an image scale 
of 0".188/pixel. For the observations during 1997 we used 
the NOT CCD#6 and for the 1998-99 observations the 
NOT CCD#7. C11601-h42 was observed in five fifiers, U, 
B, V, R and I. In the R band we also observed a neigh- 
bouring blank field, which we refer to as the background 
field. The position of this field is R.A. 16''08'"54''.0 and 
Dec. 41° 34' 00"(J2000.0). AU observations were performed 
under photometric conditions. The seeing in the images 
varies between 0".75 and l".0, depending on filter. A com- 
plete log of the observations is given in Table |^. Fig. ^ shows 
a colour image of C11601-(-42, constructed from the B, R and 
I bands. 



2.3 Data reduction 

The data were reduced using the iraf and MIDAS packages. 
Bias subtraction was made in a standard manner. Because 
of the large number of science frames in each filter, flat-fields 
were constructed directly from these images after removing 
objects. Because the /-band science frames show a fringe 
pattern we used short exposure, twilight frames for the fiat- 
field. After fiat-fielding the science /-images, a fringe frame 
was created, which was subtracted from the science images. 
The fiat-fielded images in each filter were corrected for atmo- 
spheric extinction, aligned, and finally combined. The final 
size of the field is 5'. 5 x 5'. 6 arcmin. The image s were cali- 
brated usi ng sta ndard stars from Landolt (1992) and Chris- 
tian et al. ([1983). We estimate the photometric errors in the 
zero point calibration to am ^ 0.05 in U, B, V and R, and 



0.07 in /. 



2.4 Object catalogue 



The reduced images w e re analysed using the FOC AS package 



and a minimum detection area corresponding to the seeing, 
ranging from 18 pixels in R to 30 pixels in V. The ordinary 
FOCAS procedure, including sky correction and splitting of 
multiple objects, were used. For each filter we made a cat- 
alogue including isophotal magnitude, position and area. In 
order to calculate colours of the objects we also measured 
the magnitude using a constant aperture for each object. 
Before doing this, we smoothed each image to a seeing cor- 
responding to that of the V image, which had the worst see- 
ing. The diameter of the aperture corresponds to twice the 
FWHM of the smoothed images, i.e ten pixels. A combined 
catalogue was made from the positions in the /?-catalogue 
and matching these to the positions in the other catalogues. 
A match was considered satisfactory if the distance between 
the centre positions of two objects in the different catalogues 
was less than 9 pixels. This distance was chosen to be large 
enough that each specific galaxy is found in as many colours 
as possible, but sufficiently small so that no misidentifica- 
tions occur. To check this, we visually inspected cases near 
the limit. The average offset between galaxies in the different 
colours is ;^ 2 pixels. 



2.5 Photometric limits, completeness, magnitude 
corrections and errors 

The la limiting surface brightness for the different filters 
are, in mag arcsec~^, fiu = 26.7, /ib = 27.5, /iv = 27.3, fj.R 
= 26.8 and m = 26.2. 

To estimate the magnitude correction that should be 
applied when calculating the total magnitudes from the ob- 
served isophotal magnitudes we used the same method as 
in Naslund, Fransson & Huldtgren (2000) (see also Tren- 
tham 1997). This correction accounts for the loss of light 
from the part of the galaxy profile that is outside the lim- 
iting isophote. For galaxies with exponential profiles, the 
correction is 0.04 ^ Am <^ 0.10 at /? = 25, whereas the 
correction is negligible at R < 24. For galaxies with a de 
Vaucouleurs profile, the correction is Am ~ 0.1 for galaxies 
with R < 24. We assume that only galaxies determined as 
early-types (a definition is given in Section 3.5) with /? ^ 24 
have de Vaucouleurs profiles, while fainter early-types are as- 
sumed to have exponential profiles (see discussion in Tren- 
tham 1997). We hereafter include these corrections when 
quoting magnitudes. 

To estimate the completeness, we added artificial galax- 
ies with different magnitudes, radial profiles, scale lengths 
and inclinations to the images, and used FOCAS with the 
same settings as for the real data to detect these objects. 
At R — 25, which we in the following use as our complete- 
ness limit for the photometric redshift sample, we recover 
all of the simulated galaxies. At R — 25.5 we recover ~ 97 
per cent, and at R — 26 between 69-85 per cent of the 
simulated galaxies, depending on the galaxy type. 

For bright objects the magnitude errors are dominated 
by the zero point calibration errors (given above). For faint 
objects the errors increase due to lower signal-to-noise ratio. 
At the limit, R ~ 25, we estimate the total error in R to 
be a-m ~ 0.10. 



(larvis fc Tyson 1981 



Valdes 1982 



Valdes 1993| ). We used 



a value of 3a of the sky noise as a limit to our detections, 
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Figure 1. Image of C11601+42 covering an area 3'. 7 X 3'. 7, which is ^ 50% of the total field covered. 



Table 1. Log of observations. 



Object 


Filter 


Obs. date 


Exp. time 


# of oxp. 


Seeing 


C11601+42 


U 


Apr/Jun 99 


25200s 


14 


0".89 




B 


Jun 98 


12600s 


14 


0".84 




V 


Apr 99 


16200s 


18 


1".02 




R 


Jun 97 


7200s 


8 


0".75 




I 


Jun 98/ Jun 99 


14400s 


24 


0".76 


Background 


R 


Jun 97 


7200s 


8 


0".75 



3.1 Method 

There are basically two different approaches to obtain photo- 
metric redshifts. The empirical fitting method (e.g. Connolly 
et al. 1995; Connolly et al. 1997; Wang, Bahcall & Turner 
1998; Brunner, Connolly & Szalay 1999) uses colours and 
magnitudes of objects with known spectroscopic redshifts to 
calculate a relation in form of a polynomial between these 
quantities. The observed magnitudes from galaxies with un- 
known redshifts are then inserted in this polynomial in or- 
der to calculate the photometric redshift. The alternative 
template fitting method (e.g. Puschell, Owen & Laing 1982; 
Gwyn 1995; Mobasher et al. 1996; Sawicki, Lin & Yee 1997; 
Fernandez-Soto et al. 1999; Bolzonella, Miralles & Pello 
2000; Fontana et al. 2000) uses observational or theoreti- 



cal galaxy templates of the whole range of Hubble types, 
shifted to a given redshift. A x^~^^ is made that minimizes 
the difference between the observed and the predicted fiuxes, 
magnitudes or colours. 

The advantages and disadvantages of the different 
methods are discussed i n e.g . Hogg et al. ( 1998 ), Koo ( 199!: ) 
and Bolzonella et al. (200C). A major drawback with the 
empirical fitting method is that it requires a relatively large 
set of spectra over the whole redshift range being probed. 
This is difficult, or impossible, for the faintest galaxies, and 
can easily introduce systematic errors for these. 

The weak point of the x'^~iiiethod is the template set 
used. Redshifted templates from locally observed galaxies 
are most often used to represent galaxies at earlier epochs. If 
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there are significant evolutionary changes in spectral shape 
between galaxies today and at higher redshifts, this may in- 
troduce systematic errors. A lower metallicity may have the 
same effect. Dust extinction is another source of uncertainty. 

Synthetic spectra are affected by similar effects. Internal 
dust absorption can be included by having E{B — V) and 
the extinction law as free parameters when constructing the 
template set. Effects of different metallicities can also be 
taken into account (e.g., Furusawa et al. 2000; Bolzonella et 
al. 2000), as w ell as absorption due to intergalactic H I clouds 
(Madau 1995). The latter is, however, mainly a concern for 



high-z searches, since the Lyman break is redshifted into the 
JZ-band only at ^ 2. 

Here we employ a version of the template fitting 
method. This is the most commonly used method and, de- 
spite the possible risk for systematic errors caused by an 
insufficient knowledge of high-z galaxy templates, it has 
been shown to be very successful. In the redshift range we 
are interested in, i.e., < z ^ 1.5, the e rrors are typically 
/S.Zrma/{^+z) ~ 0.05 ( Cohcu ct al. 200C ), depending ou the 



number of filters used and the accuracy of the photometry. 
The spectroscopic redshifts available for the area we have 
observed serve as a check on the accuracy of the derived 
photometric redshifts. 

For every object we minimize the expression 



i=U,B,V,R,I 



(1) 



where rrii and ai are the observed magnitudes and errors in 
the different filters, respectively. Ti{t, z) is the magnitude in 
the i-filter of template t, redshifted to z. This magnitude is 
calculated by convolving the filter curve and the quantum ef- 
ficiency of the detector with the galaxy template. The quan- 
tity ma is a constant, which fits the apparent magnitude of 
the template galaxy. 

A set of 20 different templates are used. We construct 
ten of these by interpolations between the four ob serve d 
galaxy templates given by Coleman, Wu & Weedman (1980). 
These represent E, Sbc, Scd and Im galaxies. Further, Sa 
and Sb t emplates from the K inney-Calzetti spectral atlas of 
galaxies (Kinney et al. 1996) are added to these. Finally, we 
have constructed eight templates representing E-f A galax- 
ies, using the Galaxy Is ochrone Synthesis Spectr a l Evolu- 



tion Li brary (G ISSEL) (Bruzual fc Chariot 1993; 



Bruzua] 



Chailot 1995). These templates are characterised by an 
old population (~ 8 Gyr) to which a starburst of duration 
0.25 Gyr is added. The mass of the starburst is set to be 
20 per cent of the total mass. The eight templates represent 
galaxies at epochs 0.5 Gyr to 3 Gyr after the burst. The pa- 
rameters used to construct the E-f A galaxies are cho sen t o 
be close to the preferred model used by Belloni et al. (1995), 
who compare colours from the templates with colours from a 
combination of broad and narrow band observations of E-|- A 
galaxies. 

The inclusion of E-f A galaxies are more important for 
clusters of galaxies than the field. Belloni et al. (1995) and 
Dressier et al. ( 199!: ) show that ~20 per cent of the brighter 
cluster galaxies at 0.4 ^ z <^ 0.6 have E-f A features, while 
only a few per cent of the field galaxies show these character- 
istics. Abso rptio n due to intergalactic H I clouds is treated as 
in Madau (1995). With these templates our sample include 



a fair distribution of galaxy types, including starburst and 
post-starburst galaxies. Each galaxy template is redshifted 
in steps Az = 0.002 from z = to z — 2.0. Besides the 
galaxies, we also include eight stell ar templates of M-dwarfs 
taken from Gunn & Stryker (1983). 



3.2 Accuracy of the photometric redshifts 

Dressier et al. ( 199g| ) have obtained 102 spectroscopic red- 
shifts in the field of C11601-f 42. Of these, 58 are considered 
to be cluster members, 40 belong to the field and four are 
stars. To check the accuracy of our photometric redshifts, 
we have compared our determinations with this spectro- 
scopic sample. We have excluded stars, galaxies lying outside 
our field, galaxies with spectroscopic redshifts that are not 
considered as certain (as commented in the catalogues by 
Dressier et al.) and two galaxies that appear twice in the 
MORPHS catalogues. Finally, we have excluded two galax- 
ies where the photometric redshifts together with inspection 
of the spectra (published by Dressier et al.) indicates that 
the redshifts given in the MORPHS catalogues are incorrect. 
Fig. |2| shows photometric versus spectroscopic redshifts for 
the remaining 78 comparison galaxies. The straight line in 
the figure shows the location of Zphot = Zspec- The rms devi- 
ation in redshift is ~ 0.076, and a fit to the data points 
gives Zphot = (0.98 =t 0.02) x Zspec- 

We calculate photometric redshifts for all object that 
are detected in at least three filters. If an object is not de- 
tected in all filters, we use this non-detection to set an upper 
limit to the magnitude, which is used in the redshift deter- 
mination. Of the 79 galaxies with spectroscopic redshifts we 
have magnitudes in all five filters for 63 objects, whereas 13 
objects are detected in four filters and three objects in three 
filters only. There is no difference in accuracy between the 
subsamples with detection in four and five filters. In all cases 
the [/-band magnitude is lacking. 

Liu & Green (1998) discuss the importance of the U- 
band for photometric redshift determinations. In their sam- 
ple they find that Azrms doubles if they exclude the infor- 
mation from the f/-band, compared to when they include 
all their six filters. They note that one of the main causes for 
this larger scatter is misclassification of Scd and Irr galaxies 
at z ~ 0.3 with Sbc galaxies near z ^ 0.05. At red- 
shifts z ^ 0.4, where the 4000 A break moves into the V- 
band, the importance of f/-band information is decreased. 
To confirm this we have tested this case by excluding the 
U magnitude for the 63 objects with detections in all fil- 
ters. The r.m.s. deviation is now (Jz ~ 0.11, compared to 
az ~ 0.076 with five filters. For spectroscopically deter- 
mined cluster members the accuracy is Uz ~ 0.06 with five 
filters, and ~ 0.08 with four filters. 

Because we are mainly interested in accurate redshifts 
near the cluster redshift, i.e. at 0.4 ^ z <^ 0.7, the possi- 
ble misclassification above should not infiuence the cluster 
population, and we therefore include also galaxies with only 
BVRI photometry. 

Galaxies with detection in three filters have ~ 0.05. 
The sample does, however, only consist of three galaxies and 
is therefore probably not representative. To check the three 
filter case, we instead determined redshifts for the whole 
spectroscopic sample where we only use the information in 
BRI. If we exclude the V filter, which after the U filter, is 
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Photometric redshift 

Figure 2. Photometric versus spectroscopic redshifts for 78 
galaxies in the field containing C11601+42. The cluster is located 
at z = 0.54. The deviation between the spectroscopic and the 
photometric redshifts is Azrms ~0.08. The straight line shows 
the location of Zp^ot = ^apec Spectroscopic redshifts are taken 
from Dressier et al. (1999). 



the filter reaching the smallest depth, we obtain Oz ~ 0.16. 
The increase is mainly caused by a few outliers. For spectro- 
scopically determined cluster members Oz ~ 0.13 with three 
filters. Excluding three outliers results in ~ 0.09. 

Without IR photometry the photometric redshifts in 
the interval 1 ^ z ^ 2 should be less accurate, since here 
the 4000 A break has moved long-ward of the /-band, and 
there are no other distinct features to detect. Including JHK 
photometry allows for the dete ction of the 4000 A break 
i n thi s interval. Fontana et al. (2000) and Furusawa et al. 
(2000) compares results obtained by optical four-band and 
optical+A'^Z-R 7-band photometry of galaxies in the HDF 
and find that the accuracy is considerably worse in the in- 
terval 1 < z < 2 when using optical filters only. Below 
z ^ 1 the influence of the NIR bands is, however, small. 
Since we are primarily interested in redshifts at z ~ 0.5, 
our results should not be much affected by the lack of NIR 
photometry. 



3.3 Systematic effects 

Besides the statistical uncertainty affecting the redshift 
from which the cluster membership is determined, there 
may be systematic errors in the redshift determination that 
may depend on increasing photometric errors, or even non- 
detections in some bands, at faint magnitudes, and between 
different galaxy types. 

As shown in Section 3.2, the dispersion between spectro- 
scopic and photometric redshifts increase when photometry 
is lacking in one or two filters. To examine the magnitude 
dependence of this effect, we calculate for each magnitude 
bin the mean dispersion, where we take into account the 



larger errors for galaxies that are not detected in all filters. 
In the faintest magnitude bin (24.5 < ttir < 25) we find 
that the dispersion increases by Actz < 0.006 compared to 
the over all mean. For a Gaussian distribution of the disper- 
sion, the extra fraction of galaxies lost in the faintest bin is 
< 3 per cent, clearly less than statistical errors. The rea- 
son for this is that we have chosen a magnitude limit that 
is sufficiently bright for galaxies with photometry in all five 
filters to dominate even in the faintest bin. 

Systematic effects can also be introduced because errors 
in the redshift increase at faint magnitudes due to the larger 
photometric errors. To quantify this we have made Monte 
Carlo simulations, where we create a catalog with 10^ clus- 
ter galaxies. Each galaxy is placed at the cluster redshift and 
is given magnitudes in the different bands from the spectral 
templates. To each band we add photometric errors corre- 
sponding to the magnitude dependence of the real data. We 
now determine the photometric redshift for each galaxy, and 
compare with the input galaxy redshift. In the faintest bins 
("i_R ^ 24) we find that ~ 2 per cent are outliers, i.e. have 
errors that are Az ^ 0.5. For C1160H-42 the number of 
galaxies lost due to this effect is at most a few. Secondly, 
and more important, we find that the increase in errors for 
the remaining galaxies in the faintest bins is Actz ~ 0.01, 
which translates to a loss of ~ 5 per cent. 

Using the same MC simulations we can study the error 
dependence on galaxy type. Due to the fiatter spectrum of 
the late-type galaxies, the error in the redshift is slightly 
larger compared to early-types, which have a more distinct 
4000 A break. The fraction of galaxies lost should therefore 
be higher for late-types compared to early types. All outliers 
in the faintest magnitude bins are consequently found to be 
of late-types. 

The magnitude limit of the Dressier et al. spectroscopic 
sample, which is used to calculate the errors in the photo- 
metric redshifts, is R ~ 23.3. If high-z galaxies fainter 
than this have spectral energy distributions that are not 
represented by the range of templates used, there could be 
systematic effects affecting the accuracy. The comparison of 
the spectroscopic and photometric redshifts in HDF North, 
published by Fernandes-Soto et al. (1999), does, however, 
not show any decrease in accuracy for the photometric red- 
shifts in the magnitude range 23.5 ^ R ^ 25, compared to 
the range R <, 23.5 at redshifts 0.5 < z < 2.0. This 
argues that there should be no dramatic change of the spec- 
tral energy distribution of the faint galaxies, which is not 
covered by our templates, unless faint cluster galaxies have 
a spectral shape different from faint field galaxies. 

Summarizing, we find that there is a slight systematic 
magnitude dependence of the errors in the photometric red- 
shifts, and that errors depend on galaxy type. However, the 
systematic effects introduced by these effects are unlikely to 
dominate over the statistical errors. The systematics present 
should lead to a modest underestimate of the number of faint 
galaxies, which should mostly affect late types. 

3.4 The photometric redshift catalogue 

For the objects with photometry in at least three filters we 
have constructed a photometric redshift catalogue. After ex- 
cluding 83 objects classified as stars, the catalogue consists 
of 1568 objects. 
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The catalogue contains photometric redshifts for 1140 
galaxies to i? — 25, which is 95 per cent of all galaxies to 
this li mit. Of these, 807 redshifts are based on photometry 



in five i liters, while 227 and 106 redshifts are based on four 



To estimate the number of field galaxies to our limiting 
magnitude, we have calculated photometric redshifts for the 
objec ts in the ESO Imaging Survey (EIS) ( |da Costa et al. 
199^). Using aperture photometry in U, B, V, R and I we 



and three filters, respectively. 

3.5 Determination of galaxy types 

A classification into morphological types based on visual 
(or automatic) inspection requires high resolution imaging, 
such as obtained from the HST. This is for our cluster only 
available for the central ~ 5 square arcmin. Even without 
high resolution imaging, colours can be used to estimate 
the galaxy types. This is basically what the template fitting 
method does when it compares observed magnitudes with 
magnitudes derived from a set of spectral galaxy templates. 
There is, however, not a one-to-one correspondence between 
the gal axv tvpe derived from colours and the visually deter- 



mined piorphological types. Examples are red spirals with 



little ongoing star formation ( Poggianti et al. 1999 ). Fur- 
ther, van Dokkum, Franx & Kelson (199S) find that early- 
type galaxies, especially faint SO galaxies, can have relatively 
blue colours. Also, galaxies with E-l-A type spectra can not 
be classified into any specific morphological type. Poggianti 
et al. (1999) note that galaxies in this class mainly have 
disc-like morphologies, but examples of all Hubble types are 
present. For the future discussion it is important to keep in 
mind that galaxy classification by the photometric redshift 
method is basically a spectroscopical classification, rather 
than a morphological. 

In this study we divide galaxies into two broad cate- 
gories, based on the spectral class determined by the pho- 
tometric colours of the galaxies. The division is made half 
way between the E and Sbc templates from Coleman et al. 
(1980), approximately corresponding to a rest frame colour 
B — V = 0.8. The first category consists of early-type 
galaxies with red colours, which are mostly ellipticals and 
lenticulars, but also include passive spirals with red colours. 
The second category includes late-type spirals and irregu- 
lars, with a possible inclusion of blue elliptical like sy stems . 
This division is the same as used by Ellingson et al. (2001), 
when separating red and blue galaxies. 

3.6 The redshift distribution 

The redshift distribution of the objects in the field of 
C11601-I-42 is shown in the left panel of Fig. |[ In total, 
1140 galaxies with R < 25 are included. There is a clear 
peak around the cluster location at 2 ~ 0.54, superimposed 
on a broader distribution of field galaxies with redshifts in 
the range 0.3 ^ z <, 0.7. A small secondary maximum in 
the distribution of field galaxies is apparent at 2 ~ 1.3. This 
is likely to be caused by our lack of NIR bands, since we 
above z ~ 1.2 are no longer able to detect the 4000 A break. 

A relatively large fraction of the objects in the lowest 
redshift bin is probably stars. We have already excluded 82 
objects from this bin that are best fitted by M-star spec- 
tra, or have a PSF of a point source. The magnitudes of 
the remaining objects in the bin are, however, too faint to 
allow a secure classification from an inspection of the PSF. 
The lowest bin may therefore contain a number of stars in 
addition to low redshift galaxies. 



have determined photometric redshifts for 707 objects to 
R = 25 in the 5'. 3 x 5'. 3 ElS-field including the HDF- 
South area. The redshift distribution, scaled to the area of 
C11601-(-42, of these objects is shown in the middle panel of 
Fig. 1^. Besides the peak in the distribution at 0.4 z <, 0.7, 
there is a second peak at 2: ~ 1.4, similar to the C11601-I-42 
field. Since for consistency we only include optical filters for 
the EIS data, this peak may also be caused by the lack of 
NIR photometry. 

Field-to-field variations caused by large scale structures 
introduce an uncertainty when correcting for field galaxies. 
To check if the EIS field has a number density that is com- 
parable to other fields, we use estimates from Small et al. 
([l997|), w ho find that the Medium Deep Survey ( [Sriffiths et 
al. 1994) has a field galaxy density of ~ 11.2 galaxies per 
"arcmhF to I = 22.6. From the EIS field we find ~ 11.0 
galaxies per arcmin^ to the same magnitude limit suggest- 
ing that the background counts in the EIS field are normal. 
Deviations may, however, increase at fainter magnitudes. To 
quantify the field-to-field variations, we calculate the vari- 
ance in the background counts from the angular correlation 
function (Peebles 1975; 1980). At J? > 24 we use a correla- 



tion amplitude from Connolly et al (1998), while for brighter 
magnitudes we calculate the correla tion a mplitude from a fit 
to the data in Woods & Fahlman (1997). The la variation 
is found to be ± 13 per cent at ~ 25 for our field size. 



3.7 Cluster membership 

The limited accuracy of the photometric redshifts makes 
the determination of cluster membership somewhat compli- 
cated. First, a redshift range defining the cluster has to be 
determined. Then, the contamination of field galaxies within 
this redshift range should be estimated. 



For the spectroscopic sample. Dressier et al. (1999) use 
a redshift range 0.5100 < z < 0.5474 to define cluster 
membership. In the photometric sample the higher disper- 
sion in redshift has to be taken into account, and a broader 
interval must be used. The redshift interval has to be large 
enough so that as many of the cluster galaxies as possible are 
retained, but at the same time keeping the number of field 
galaxies in the cluster redshift range to a minimum. Here, 
we take a range around the mean redshift of the cluster cor- 
responding to 1.5az, corresponding to zci= 0.43-0.65. The 
interval is shown in the right panel of Fig. ^ and contains 
418 galaxies. Correcting for the few percent of galaxies with 
R < 25 that lack photometric redshifts results in a total 
of 432 galaxies. After renormalising for the different size of 
the EIS field, we find that 185 field galaxies are expected 
within the above redshift range. The total number of cluster 
galaxies after subtracting the background within the cluster 
redshift range is therefore 247. 

The 1.5o" cut in the redshift range means that we delib- 
erately loose some cluster galaxies. Increasing the redshift 
range would, however, increase the background contamina- 
tion, and therefore increase the uncertainty. In particular, 
the probability that a specific galaxy is a cluster member 
would decrease. For a Gaussian distribution of the redshift 
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Figure 3. Left panel: Distribution of photometric redshifts for 1140 objects in the field containing C11601+42 
to magnitude R < 25. Middle panel: Same for the ESO Imaging Survey, corrected to the same field size 
as the C11601+42 field. Right panel: Residuals after subtracting the background field from the C11601+42 
field. The redshift range chosen for the cluster analysis is marked by the horizontal bar. 



errors we estimate that ~ 20 per cent of the cluster popu- 
lation is outside the redshift range. Here, we have assumed 
that the r.m.s. error in the catalogue is higher than in the 
comparison with the spectroscopic sample by Act — 0.01, 
since a slightly larger fraction of galaxies in the catalogue 
do not have photometry in all five filters, as compared to 
the spectroscopic sample. When appropriate, we correct for 
this loss of galaxies in the analysis. 

The difference between using photometric redshifts and 
an ordinary background subtraction is that we only have to 
do a subtraction for a limited redshift range. This decreases 
the errors compared to using the full range of redshifts. A 
simple estimate, based on Poissonian statistics, illustrates 
this. Of the total number of ~ 1200 galaxies we estimate 
that 247 are cluster galaxies and the remaining ~ 950 are 
field galaxies (for simplicity we do here not correct for galax- 
ies lost due to the 1.5a cut). With a standard subtraction 
method, the r.m.s. error in the number of cluster members 
is equal to the r.m.s. error in the subtracted background, 
±\/950. This makes the number of cluster galaxies 247 ±31. 
When only subtracting the ~ 185 field galaxies within the 
cluster redshift range the r.m.s. error is reduced to ± y/TEE, 
and the number of cluster galaxies becomes 247 ± 14. The 
r.m.s. error is therefore reduced by more than a factor two. 

Most important is that the probability that a specific 
galaxy is a cluster member increases from 247/1200 ~ 21 
per cent without redshift information to 247/432 ~ 58 per 
cent with redshift information. With a field contamination as 
high as ~ 80 per cent, it is impossible to obtain useful infor- 
mation about e.g. the radial distribution within the cluster 
in the standard subtraction method. With the limited sub- 
traction necessary when using the photometric method, we 
can, however, derive meaningful information about the in- 
ternal properties of the cluster. In combination with the rich 
colour information of the individual galaxies, this method al- 
lows more complex investigations of the cluster properties. 
A further advantage is that the field-to-field variations are 
minimized. 



4 RESULTS 



4.1 Cluster properties 



Abell (L958) defines Richness classes by the number of 
galaxies in the cluster with magnitudes between ma and (ma 
-|- 2), where ms is the third brightest galaxy in the cluster. 
Here we find ms ~ 20.1 in the i?-band, corresponding to 
Mb — — 22.5. We find 135 galaxies in this interval. Af- 
ter subtracting for background galaxies within the redshift 
range of the cluster (see below), we estimate the resulting 
number of cluster galaxies to ~ 70. This makes C11601-I-42 
a cluster of Abell class 1-2, i.e. an intermediate rich clus- 
ter. This is important to keep in mind when comparing with 
other results, since most investigations concern rich clusters. 

Clusters can also be characterised by the degree of cen- 
tral concentration. Butcher & Oemler (1978b) define a con- 
centration index, C = log(i?6o/^2o), where R„ is the radius 
of a circle containing n per cent of the projected galaxies 
in the cluster with My < — 20 (using a cosmology with 
qo = 0.1 and Ho = 50 km s~^ Mpc'^). B084 define 
clusters with C ~ 0.3 as open or irregular. Clusters with 
C <; 0.4 are normally referred to as compact or concentrated. 

We estimate that C11601-I-42 has C ~ 0.40 ± 0.04. 
Because our field does not contain the whole cluster, we 
have extrapolated the radial profile according to an isother- 
mal sphere outside of our image. This mak es the value of C 
somewhat uncertain. Dressier et al. (1997) find C = 0.34 
for C11601-I-42. The difference in C, if significant, is proba- 
bly due to different methods used for calculating the back- 
ground contamination. Of the clusters in the MORPHS sam- 
ple, C1160H-42 belongs to the group with relatively low con- 
centration. 



4.2 The cluster luminosity function 

4-2.1 The photometric redshift method 

To calculate the LF we divide the 418 galaxies within the 
cluster redshift range into magnitude bins of Am = 0.5 to 
R — 25. We smooth each galaxy count with a Gaussian 
using a a corresponding to the error in magnitude. The effect 
of this is small since the errors are typically much less than 
the bin size. It does, however, result in bins having fractions 
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Figure 4. Luminosity function in the ij-band for C11601+42. 
Error-bars include Poissonian uncertainties and field-to-field vari- 



of galaxies, which can be noted in the brightest bins. In the 
faintest magnitude bins we correct for the fraction of the 
galaxies that lack photometric redshifts by multiplying the 
number of galaxies with the inverse of the completeness of 
the photometric catalogue. After correcting for field galaxies 
in the redshift range of the cluster for each magnitude bin, 
we obtain the LF shown in Fig. Q The error-bars in the 
figure represent la, including Poissonian uncertainties and 
field-to- field variance. 

The solid line in the f igure shows a fi t to the LF by the 
usual Schechter function dSchechter I976| ), 



$(M) oc e 



10 



-0.4(a+l)A/ 



(2) 



where M* is the characteristic magnitude, representing the 
turnoff at the bright end of the LF profile, and a is the 
slope beyond the up-turn at the faint end of the LF. We 
find a = — 1.42 ± 0.12, and a characteristic appar- 
ent i?-band magnitude m|j = 20.8 ± 0.6. The errors 
in the Schechter function parameters are calculated from 
Monte Carlo simulations, where we add, or subtract, a ran- 
dom number of galaxies to each magnitude bin according to 
the uncertainties. The quoted errors are the la deviations in 
the parameters from one hundred Monte Carlo realisations. 
The LF for the background field at the cluster redshift has 
= 20.6 ± 0.2, and a = -0.90 ± 0.08. This is similar 
to the slope for the tot al LF for the field, a ~ — 1, found 
by e.g . Loveday et al. (1992) and Marzke, Huchra & Geller 

(Fggi). 



The determination of M* is based primarily on the 
brightest galaxies. Here statistics are poor, and the error 
in M* significant. For example, excluding the two brightest 
bins, with only three galaxies, increases M* by ~ 0.6 mag. 
There is also a strong coupling between M* and a, making 
the faint-end slope dependent on the derived value of the 
turnoff magnitude. In Section 5.1 we show that a Gaussian 



plus Schechter function fit for the cluster LF gives a better 
representation of the underlying population. To avoid the 
coupling between M* and a we also calculate the faint-end 
slope by fitting the five faintest magnitude bins, i.e. ~ 2.5 
mag brighter than the limiting magnitude, to a straight line 
(as proposed by Trentham 1998). The line has the form 



(3) 



where the index / denotes that only the faint magnitude bins 
are used. For C11601-I-42 this yields a; = - 1.38 ± 0.07. 
For the background field we find a/ = — 1.17 ± 0.10. 

In order to compare our results to the slope for other 
more nearby clusters, we must add K-corrections to calcu- 
late the absolute rest frame magnitudes in different bands. A 
major advantage of the photometric method is that individ- 
ual K-corrections for each galaxy are known, based on the 
best-fitting template. In the standard one-filter subtraction 
method, a single value for the K-correction usually has to 
be assumed, due to the difficulty of determining the galaxy 
type of the faint objects from single band observations (e.g. 
Driver et al. 1998a; Naslund et al. 2000). 

In the left panel of Fig. ^ we plot the resulting rest frame 
_B-band LF. For a constant K-correction this should be sim- 
ilar to the observed 7?-band LF. Because of the variation in 
the individual K-corrections of the different galaxy types 
there is, however, a small redistribution of the galaxies be- 
tween the bins, leading to a slight steepening of the faint-end 
slope. The turnoff magnitudes, faint-end slopes and abso- 
lute magnitude limits for the rest frame LFs in all bands are 
listed in Table ^ The absolute magnitude limits for the dif- 
ferent filters approximately correspond to R = 25, and are 
calculated from the average of the individual K-corrections 
for the galaxies. 

4-2.2 The cluster LF for different populations 

To investigate the LF of the different galaxy types, we divide 
the cluster galaxies into early-type and late-type subsam- 
ples, as described in Section 3.5. It is important to note that 
this refers to a spectroscopic type, not morphological. Field 
galaxies within the cluster redshift range are subtracted with 
a distribution of galaxy types determined from the EIS data. 
The late-type fractions in the cluster to different limiting 
magnitudes are listed in Table Early-types dominate at 
bright magnitudes, whereas late-types dominate at magni- 
tudes somewhat fainter than R ~ 23. 

The distinction between a blue population that domi- 
nates at faint magnitudes, and a population of bright red 
early-type galaxies is illustrated in Fig. ^ where we plot 
the total rest frame _B-band LF, as well as the LFs for 
the two populations separately. The faint-end slope and 
characteristic magnitude for the late-type subsample are 
Q = - 1.68 ± 0.11 and M*b = - 21.4 ± 0.6. For 
the early-type subsample we find a — — 0.41 ± 0.20 and 
Mg = - 21.3 ± 0.4. The low statistics at Mb <, - 20.5, 
of the late-type LF makes this part, as well as Mg, un- 
certain. The straight line fit for Ms ^ — 20.1 results in 
Q/ =-1.80± 0.14, and a/ = 0.17 ± 0.42, for the late- 
type and early-type subsamples, respectively. 

It is clear that the faint-end slope of the late-type LF 
is steep, and that this blue population also dominates the 
faint end of the common LF. The red population shows an 
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Figure 5. Luminosity functions in the rest frame _B— band for C11601+42. The left panel shows the total 
LF, whereas the middle and right panel shows the LFs where the galaxies are divided into early- types and 
late-types. The total and late-type LFs are fitted to a Schechter function. The solid line in the middle panel 
shows a Gaussian fit, while the dashed line is based on the Schechter function. Note the better fit of the 
Gaussian. 



Table 2. Parameters for the LF in the rest frame for different filter bands. M* and a are the characteristic 
magnitude and faint-end slope of the Schechter function, while cm* and Oa are the errors in these parameters. 
af IS the slope of the straight line fit to the five faintest magnitude bins, and CTq^ is the error in this quantity. 



Filter 


Magnitude limit 


M* 


am- 


a 


(Tq 


aj 




U 


-17.9 


-21.8 


0.9 


-1.46 


0.25 


-1.49 


0.12 


B 


-17.6 


-22.3 


0.5 


-1.37 


0.07 


-1.41 


0.09 


V 


-18.4 


-22.5 


0.8 


-1.40 


0.14 


-1.44 


0.09 


R 


-18.9 


-24.0 


0.9 


-1.48 


0.11 


-1.49 


0.09 


I 


-19.5 


-24.8 


0.9 


-1.41 


0.11 


-1.45 


0.10 



Table 3. The fraction of the total number of galaxies that are 
late-type spirals or irregulars at different limiting magnitudes in 
the _R— band. 



Magnitude, 

Late-type fraction 



< 23 23-24 24-25 



0.41 



0.89 



0.96 



opposite behaviour at faint magnitudes, and is better fitted 
by a Gaussian. The reduced chi-square of tiie Gaussian fit to 
the early- type LF is x^/i^ ~ 0.54, while tiie Sciiechter func- 
tion fit yields "x^ /v ~ 1.0, clearly showing that a Gaussian 
gives a better representation of the LF. 



4-2.3 The standard background subtraction method 

As a comparison, we calculate the LF using the standard 
background subtraction metliod in a single filter. As back- 
ground we use our _R-band image of a blank field, and sub- 
tract this from the cluster image, correcting for the slightly 
different field sizes. With this method we reach R = 26, 
i.e., about one magnitude fainter than for the photometric 
redshifts. For the two faintest magnitude bins we have cor- 
rected for incompleteness as described in Section 2.5. This 
results in 552 ± 42 cluster galaxies to i? = 26, while the 
number of galaxies to i? = 25 is 373 ± 27, where the 



errors are statistical only. Including field-to-field variance in 
the background counts results in 373 ± 120 at J? = 25. 
To compare this with the number found by the photometric 
method, we correct the latter for the galaxies lost due to the 
1.5a cut in redshift range, and include field-to-field variance. 
This results in 309 ± 31 at J? = 25, which is within the 
errors found by the subtraction method. 

The LF to 7? = 26 using the subtraction method 
is shown in Fig. ^. A Schecter fit yields a faint-end slope 
Q ~ — 1.38 ± 0.12, and the turnoff magnitude is 
m|j = 20.6 ± 0.3, close to the values derived using 
photometric redshifts. The faint end fit, as defined by Eq. 
(3), for both Rum ~ 25 and Rum = 26 results in 

~ — 1.44 ± 0.11 which shows that there is no dramatic 
change at the faint end when going slightly deeper in mag- 
nitude. Comparing the error-bars in Figure ^ and Figure g 
clearly illustrates the gain with the photometric method, es- 
pecially at faint magnitudes. The larger error-bars in Figure 
W around R ~ 21 are due to the smaller statistics in these 
bins for the photometric method, probably caused by the 
use of different background fields in the two methods. 

The result also indicates that the subtraction method 
works reasonably well to redshifts z ~ 0.55, even though 
the number of cluster galaxies vary somewhat between the 
method s. This result may seem contrary to Driver et al. 
(1998a), who from simulations find that the method is un- 
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Figure 6. Luminosity function in the i?-band for CI1601+42 cal- 
culated by the standard method of subtracting a background field 
from the cluster image. Error-bars include Poissonian uncertain- 
ties and field-to-field variance. 



reliable for clusters ai z > 0.3. Their conclusions are, how- 
ever, valid for a specific observational set-up. We reach a 
completeness limit ~ 1.5 mag fainter than the limit used 
by Driver et al., which should increase the reliability of the 
subtraction method to a higher redshift. 



4-2.4 Gravitational lensing and crowding in the 
subtraction method 

A possible source of uncertainty in the background subtrac- 
tion method comes from gravitational lensing. Background 
galaxies in the cluster image can have their magnitudes am- 
plified by the cluster, leading to an excess of background 
galaxies in the cluster image, compared to the cluster free 
image used for subtraction. For the rich cluster A963, which 
show signs of arcs, Driver et al. (1994) estimate that the 
increase in background counts is ~ 10 per cen t. W hen in- 
vestigating the Coma cluster, Bernstein et al. (1995) argue 
that this increase will be cancelled by the reduction of the 
solid angle being viewed behind the lensing cluster, mak- 
ing no correction necessary. Considering this, together with 
the fact this C11601-I-42 is both intermediate in richness and 
concentration and has no detected arcs, we find that grav- 
itational lensing should not affect our results. Effects from 
gravitational lensing should not enter into the photometric 
method, since we do not subtract galaxies behind the clus- 
ter, but instead reject them using the redshift information. 

Another factor to account for in the subtraction method 
is crowding, i.e. overlapping object. Especially faint objects, 
with small areas, can be hidden behind brighter galaxies and 
stars. Since there are in total more objects in the cluster 
image than in the background image, this can lead to a 
higher fractional loss of faint galaxies in the cluster image, 
resulting in a LF with too few faint objects. Using similar 



arguments as in Naslund et al. ( ^000[ ), we estimate the loss 
by adding the total area covered by objects in the cluster 
image and the background image, respectively. We find that 
the correction in the number of faint galaxies in the cluster 
LF is less than two per cent. Since this correction factor 
is much less than statistical errors, as well as field-to-field 
variations, we do not include any corrections due to crowding 
in our results. 

For the photometric method, only foreground galaxies 
contribute to the crowding. The effect should be of equal 
importance in the cluster field and the field used for sub- 
traction, and therefore cancel. 



4.3 Radial distribution 

Using the cluster galaxies selected by photometric redshifts, 
we plot the radial distribution of the cluster members with 
i? < 25 in Fig. ^. The left panel shows the projected surface 
density, and the right panel shows the fraction of late-type 
and early-type galaxies as a function of radius. 

The most interesting result is that the cluster outside 
of ~ 0.7 Mpc shows a roughly constant fraction of late- type 
galaxies, ~ 70 per cent of the total. Inside this radius the 
fraction of early-type galaxies increases rapidly to ~ 80 per 
cent at the centre. Although, this effect is present also in low 
redshift clusters (e.g. Whitmore, Gilmore & Jones 1993), our 
result shows that this population gradient is set up already 
at z ~ 0.54. In fact, the gradient is even steeper than in 
local clusters. 

In Fig. ^ we plot the radial distribution for bright 
{R < 23) and faint (23 < R < 25) galaxies separately. 
While the bright galaxies show a steep cluster-centric dis- 
tribution, the faint galaxies have a slower increase, which 
flattens at the centre, consistent with a homogeneous distri- 
bution of faint galaxies in the cluster. The figure also shows 
the radial distribution separately for bright and faint late- 
type galaxies. We discuss the radial distribution further in 
Section 5.2. 



4.4 The Butcher— Oemler effect 

B084 showed that the fraction of blue galaxies, fs, in 
cores of rich compact clusters increases dramatically from 
/s ~ 0.03 ± 0.01 for clusters with z ^ 0.1, to /s ~ 0.25 
at 2 ~ 0.5. 

To determine /s in C11601-I-42 we adopt the definition 
in B084. First, galaxies brighter than Mv = — 20 within 
a circular radius, -R30, containing 30 per cent of the cluster 
galaxies are selected. From this sample fn is defined as the 
fraction of these galaxies with rest frame colours B — V 
at least 0.2 mag bluer than the linear locus of the colour- 
magnitude (CM) relation. In order to compare our results 
with those of B084 we use in this section their adopted 
cosmology, i.e Ho = 50 km s~^ Mpc~^ and go = 0.1. 

We find _R3o ~ 1.2 arcmin, or ^ 0.6 Mpc for the 
cosmology above. This is somewhat smaller than the t ypica l 
radius of R30 ~ 1 Mpc given by Kodama & Bower (2001) 



for a sample of rich clusters at 0.2 <^ z <^ 0.4, but similar 
t o the radius R30 ~ 0.5 Mpc used by van Dokkum et al. 
( ^OOOl ) for the cluster MS1054-03 at z = 0.83. Our value 
of i?3o is also similar in angular size, ~ 1 arcmin, to that 
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Figure 7. The left panel shows the radial distribution of the cluster galaxies in C11601+42. Background 
subtraction within the cluster redshift range has been applied. The right panel shows the fraction of early- 
type and late-type galaxies in the cluster as a function of radius. In absolute magnitudes _R < 25 corresponds 
to Mg < — 17.6 in the rest frame. 
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Figure 8. The left panel shows the radial distribution of bright (R < 23, Mb ^ — 19.6) and faint 
(23 < R < 25, -19.6 ^ Mb ^ - 17.6) galaxies in C11601+42. The right panel shows the radial 
distribution of bright and faint late- type galaxies. 



used by B084 for their highest redshift clusters. The blue 
fraction we find is /s — 0.15 ± 0.06. 

This value for C11601+42 is significantly higher than 
the local value /s ~ 0.03, but lower that the average value 
/b ~ 0.27 &t z = 0.54, taken from the fit presented by 
B084. Note, however, that B084 did not use K-corrections 



that de pend on the spectral type of the galaxies ( Kodama & 
Bower !!001). Including this, /a decreases by ~ 0.04, and 
the K-corrected average BO value at 2: = 0.54 would there- 



fore be /s ~ 0.23. Even after this correction, we find that a 
significant difference between the blue fraction in C11601-I-42 
and the B084 sample remains. 

Finally, a possible concern here are systematic ef- 
fects, which can discriminate against selecting blue late-type 
galaxies, and therefore lead to an underestimate of /b. From 
the MC simulations described in Section 3.3, as well as com- 
paring with available spectroscopic redshifts, we find that 
systematic effect can cause a relatively higher loss of blue 
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galaxies when selecting cluster galaxies by using a fixed red- 
shift interval. This may lead to an underestimate of the blue 
fraction by A/s ~ 0.01. The small effect here is due to the 
relatively bright limiting magnitude used to select the BO- 
galaxies, i.e. Mv < — 20 (corresponding to niR <^ 23.5). 
Here systematic effects are not expected to affect the pho- 
tometric redshifts more than marginally. However, when we 
in Section 5.3 determine the blue fraction to fainter magni- 
tudes, the underestimate may reach Afs ~ 0.05. 



4.5 Comparison with Hubble types from HST 

C11601+42 was observed with WFPC-2 on HST during Cy- 
cle 4 with the F702 W fiher. The total exposure was 16800 
s (Small et al. 1997). The image covers the central 5 square 



arcmin of the cluster, and detailed morphological informa- 
tion are presented by Small et al. for all 145 objects with 
i?F702w < 23.5. Five of these are classified as stars, whereas 
no classification was possible for a further six objects. Of the 
five objects classified as stars we get the same identification 
for three, while the other two are found in the lowest redshift 
bin, consistent with these also being stars. Of the six objects 



Of 



withoul HST classification we identify three as stars, one as 
possiblf star and two as galaxies. This leaves 134 galaxies 
with classified morphology. We are able to determine pho- 
tometric redshifts for 129 of these objects. 



the 68 gfllRTies classified as parly-type galayies by 



the photometric redshift method, ~ 62 per cent have mor- 
phologies classified as E to Sa, and the rest are Sab to Irr. 
Conversely, of the 61 galaxies photometrically classified as 
late-type galaxies, ~ 72 per cent are morphologically clas- 
sified as Sab to Irr. Even though there may be some uncer- 
tainty in the visual classification into Hubble types (Small 
et al. set Rf702w < 23.5 as the limit where no more than 
20 per cent of the galaxies differ by more than one Hub- 
ble type between two different classifiers), the main reason 
for this relatively low correspondence is the lack of a direct 
relation between spectroscopic and morphological types, as 
also show n by e.g.. Dressier et al. (1999) and Poggianti et 

al. (Fgggh. 



This can be seen in Fig. 9, which shows the rest frame 
B — V colour for the galaxies within the redshift range of 
the cluster, divided into two broad morphological classes. 
The lower panel shows the colours of the galaxies with early 
Hubble types (E-Sa), and the upper panel later types, as de- 
termined from the HST data by Small et al. Most early-type 
galaxies are red and have a small dispersion in colour. There 
is, however, a tail of galaxies with relatively blue colours. 
The distribution of later types have a broader distribution, 
but are on average ~ 0.2 mag bluer than the earlier types. 
The figure clearly shows that there is relatively large popula- 
tion of objects, morphologically classified as late-type galax- 
ies, that have red colours. 



5 DISCUSSION 

5.1 Evolution of the cluster luminosity function 

Luminosity functions for galaxies fainter than Mb ~ — 18 
have previously been determined for clusters to z ^ 0.3. 
Comparing the LFs for Coma at 2 — 0.023 (Trentham 
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Rest frame B — V colour distribution of early-type 
-Sa) and late-type galaxies (Sab— Im) brighter than 
23.5 in C11601-I-42. Morphologies are as determined 
al. (1997). 



1998), A2554 a.t z = 0.106 (Srn ith et al. 1997), A963 at 
z = 0.206 ( [Driver et al. 1994| ) and MS2255. 7-^2039 at 
z = 0.288 (Naslund et al. 2000) reveals that both nearby 
and distant clusters show a steep faint-end slope (q <^ — 
1.4), indicating no qualitative difference with redshift. Our 
results for C11601-I-42 show that this conclusion is true up 
to z ~ 0.5. 

For the highest redshift cluster with a deep LF, 
MS2255.7-h2039 a.t z = 0.288, Naslund et al. ( |200C| ) de- 
rive M|j = — 24.2 (in our cosmology) and a ~ — 1.43, 
when fitting their data to a single Schechter function. This 
agrees well with our results in Table |^ 

The increased fraction of faint galaxies in the outer part 
of C11601-I-42, shown in the left panel of Fig. |, has the 
consequence that the faint-end slope of the LF should de- 
pend on radius. Consequently, we find that the faint-end 
slope for galaxies within a radius of 80 arcmin (0.7 Mpc) is 
Uf ~ - 1.27, while between 80-160 arcmin (0.7-1.4 Mpc) 
a/ ~ — 1.52. When comparing the LF's between clusters, 
it is therefore important to remember that the slope of the 
LF in the core region can be much shallower than for the 
cluster as a whole. 

A steepening of the slope with r adius is consistent with 
studies of local clusters. Adami et al. (2000) e.g., find a faint- 
end slope a ~ — 1.0 for the c ore r egion of Coma, while 
for the whole cluster Trentham (1998) finds a considerably 
steeper faint-end slope, a ~ — 1.7. Also, five of seven clus- 
ters at z ~ 0.15 examined by Driver, Couch & Phillipps 



( 19981 ) show this trend. The faint-end slopes for these seven 



clusters vary between -0.9 and -1.4, which is more shallow 
than the slopes for the clusters quoted in the beginning of 
this section. The flatter slope compared to C11601-I-42 may 
be understood as a result of the effect above, since the clus- 
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Figure 10. B-band LFs for C11601+42 and Coma (from Tren- 
tham 1998). Note the different sliapes of tlie LP at intermediate 
magnitudes -20.5 ^ Mb ^ - 18.0. At the faint end of C11601+42 
we have added two data points derived from the single-filter sub- 
traction method (shown as open circles). Coma is normalised to 
C11601+42 at Mb ~ - 20. 



ters in the sample of Driver et al. only cover a radius half 
the linear size of our cluster. 

It is interesting to compare the rest frame B-band 
LFs of C11601-I-42 and the Coma cluster (Trentham 1998), 
shown in Fig. Coma is here normalised to C11601-I-42 at 
Mb ~ — 20. Note that we use similar methods as Tren- 
tham to correct the magnitudes at faint levels (see Section 
2.5). Also, both investigations use a 3(t detection and similar 
surface brightness limits. We have for C11601-I-42 added two 
data points derived from the one colour subtraction method 
(open circles). These points are in agreement with a con- 
tinuation of the steep slope found at intermediate magni- 
tudes. The figure shows that at intermediate magnitudes, 
—20.5 <, Mb ^ — 18.0, Coma has a Gaussian LF peaking 
at Mb ~ — 19.5, before the faint galaxies start to domi- 
nate at Mb ^ — 18.0, and the slope of the LF turns steep. 
However, for C11601-I-42 the steep part sets in already at 
Mb ~ — 20, and no plateau is present, as in Coma. 

If the apparent lack of a 'plateau' in C11601-I-42 is a 
general property at this redshift, it could indicate a general 
brightening of the faintest cluster members, because of more 
active star formation, with increasing redshift. This bright- 
ening would shift the faint, steep end towards brighter mag- 
nitudes at higher redshift, and therefore decrease the extent 
of the 'plateau'-range, while preserving the slope. This is 
consistent with Coma having a/ ~ — 1.5 at Mb ^ — 17, 
and C11601-h42 a/ ~ - 1.4 at Mb k, - 19. 

The 'plateau' between the turnoff at bright magnitudes, 
and the up-turn at faint, is present in both MS2255 at 
z = 0.288 and A963 a.t z = 0.206. This indicates that 
the star formation activity of the faint galaxies indicated by 



C11601-I-42 at z = 0.54, is not as important in the two 
clusters at lower redshift. 

There is some controversy regarding the colour depen- 
dence of the faint-end slope of the LF. De Propris et al. 
( [1995| ) determine the LF for four Abell clusters in both B 
and /. Although the slopes are very steep (a < — 2), 
there is no difference between the two colours. These results 
are in contrast to the results of Wilson et al. (1997), who 
determine the LF in B and / for A665 and A1689, both 
at z — 0.18, and find a ~ — 2 for the y-band and 
a ~ — 1.2 for the /-band. Comparing the slope of the LF's 
of the different colours for C11601-I-42 (Table ||), there is no 
significant trend when going from red to blue filters. This is 
consistent with the dominance of the late-type population 
at faint magnitudes. If a majority of the galaxies have simi- 
lar colours, i.e. K-corrections, then one expects a horizontal 
shift of the faint-end LF between the different bands, while 
the steepness remains unchanged. 

To investigate if there is a trend at even shorter 
wavelengths in C11601-I-42, we calculate a/ at rest frame 
UV{2850 A) and UV{2350 A), corresponding to the ob- 
served B- and [/-bands, and find a/2850 = — 1.71 ± 0.20, 
and Of 2350 = - 2.05 ± 0.20, for Muv (2850) ^ -19.0 and 
Muv(2350) Z -19.1, respectively. 

While the optical bands do not show any trend in the 
faint-end slope, the steepening of the slope in the UV-hands 
suggests that the faintest cluster galaxies have an increasing 
fraction of young, UV bright stars. 

In a recent paper Valotto, Moore & Lambas (2001) dis- 



cuss systematic eflects when determining the LF for clus- 
ters selected in two dimensions, such as the Abell and 
EDCC (Edinburgh-Durham Cluster Catalogue, Lumsden et 
al. 1992) catalogues. From N-body simulations Valotto et. 
al. find a bias towards selecting clusters which are contam- 
inated by background clusters, or groups, projected on to 
the cluster. The over-density in the background, along the 
cluster line-of-sight, will result in an over-abundance of faint 
galaxies, which remains after the subtraction of a cluster-free 
background field. For a composite LF, calculated from sim- 
ulations where the clusters are selected in two dimensions, 
similar to the Abell clusters, Valotto et al. find that the selec- 
tion process therefore can introduce an artificial steepening 
of the LF. 

By selecting cluster members with photometric red- 
shifts, the background contamination should be minimized 
in our sample. The steep faint-end slope that we find should 
therefore not be much affected by the selection effects de- 
scribed by Valotto et al. Furthermore, the similar slope we 
derive from the background subtraction method indicates 
that projection effects are not important for C11601-I-42. 



5.1.1 Evolution of the LF for different cluster populations 

From deep studies of both near and intermediate redshift 
clusters, it is clear that a single Schechter function gives 
a rather poor representation of the LF (e.g. Driver et al. 
1994; Wilson et al. 1997; Naslun d et al. 2000). For the 



Virgo cluster Binggeli et al. ( 1988 ) find that ellipticals with 
Afs <, — 18, as well as the different types of spiral galaxies, 
are better represented by Gaussian LFs. The mean magni- 
tude is for the spirals shifted towards fainter magnitudes for 
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the later types. The sum of all types of spirals can also be 
well represented by a Gaussian LF. 

At magnitudes Mb ^ — 18 dwarf ellipticals (dE) and 
irregulars dominate. The up-turn of the LF at the faintest 
magnitudes is mostly due to the dE population, where a 
Schechter function is a better representation of the LF. The 
total LF should therefore rise as a Gaussian at bright mag- 
nitudes and then flatten to a 'plateau' at intermediate mag- 
nitudes, before it again turns steep at the faint end. 

The conclusions above may seem to contradict other 
studies showing a flat LF. Garilli, Maccagni & Andreon 
(1999) flnd a ~ — 0.9 for a composite LF calculate d from 
65 clusters at 0.05 < z < 0.25, while PaohUo et al. ( |200l| ) 
flnd a ~ -1.1 for 39 clusters at 0.058 < z < 0.28, when 
fltting their data to a single Schechter function. Paolillo et 
al. do, however, point out that the LF below M ~ — 19.5 in 
these clusters, which probes the dwarf population, is system- 
atically above the Schechter function fit, indicating a upturn 
of the faint-end slope at these magnitudes. The flat slopes 
found in these investigations are therefore not in conflict 
with the steeper slopes found in deeper studies. 

Above we argued that the lack of a plateau in the LF 
of C11601-(-42, as compared to clusters at lower redshift, 
can be explained by a uniform shift in magnitude for the 
fainter galaxies. Because these galaxies are predominantly 
blue (Fig. |), the proposed shift can be interpreted as a 
result of a higher star formation level in this population 
of galaxies compared to the bright population. For cluster 
galaxies at z = 0.24, which have had their star formation 
truncated, Small et al. (1998) show that late types galaxies 
both fade more and become redder than early type galaxies. 
This is consistent with the shift of the late type population 
in C11601-I-42. The total reddinging and fading between z = 
0.24 and z = is, however, fairly modest in these models. A 
more substantial fading is expected if dynamical processes 
are involved. 

Wilson et al. ( 1997| ) argue that dwarf irregulars in clus- 
ters should fade rapidly, either due to a trunction process, 
such as ram pressure stripping, or by disruption due to tidal 
interactions. They estimate the fading since z ~ 0.2 to ^ 3 
mag. Such a fading would make part of the progenitors of 
todays dwarfs with Mb i^, -15 detectable in C11601-I-42, and 
could therefore contribute to the faint-end slope. Further 



Kodama fc Bower 200l| ). A bluening with redshift, i.e. an 
increased star formation rate, of the infalling held popu- 
lation can also contribute to the observed effect. Diaferio 



simulations of 'galaxy harassment' (Moore, Lake & Katz 
1996) Slow that high redshift {z ~ 0.4) cluster galaxies of 
type Sc or later can be turned into low surface brightness 
spheroidal systems in todays local clusters. The fading of 
these 'harassed' galaxies should be ~ 2 mag. The dwarf el- 
lipticals present at z ~ 0.5 should, on the other hand, only 
fade marginally by 0.5 mag. Whatever the mechanism, it 
seems that at least a fraction of the galaxies contributing to 
the steep slope at the faint end in local clusters have faded 
considerably compared to earlier epochs. The gradual dis- 
appearance of the plateau at high redshift could therefore 
be explained as a fading of the faint, late- type population, 
relative to the bright population. 

This scenario requires that the number of recently ac- 
creted galaxies is higher at z ~ 0.5 than today, i.e. that 
there is an increase in the infall rate with redshift. As shown 
by the Press-Schechter formalism, N-body si mulations, as 
well as observations, this is indeed the case (IBower 1991 



et al. (2001) show that there is a significant increase in 
star formation rate in field galaxies over the redshift range 
0.18 < z < 0.55. Therefore, the brightening is likely to be 
due to a combination of the above two effects. The bright- 
ening of the blue cluster population also relates to the BO 
effect, which shows an increase of star formation in a large 
fraction of cluster galaxies at higher redshifts. We return to 
this in Section 5.3. 

We now focus on the e volution of the bright cluster pop- 
ulation. Small et al. (1997) calculate composite LFs for ten 
clusters in the MORPHS sample divided into three redshift 
bins centred on z — 0.38, z — 0.43 and z — 0.54. (The last 
bin includes C11601-I-42.) The galaxies are here divided into 
ellipticals and spirals (Sa and later) , based on visual classifl- 
cation into morphological types. The limiting magnitude for 
this sample is, however, too bright to detect the dwarf popu- 
lation, or for determining the faint-end slope. Instead, when 
calculating the turnoff magnitude. Small et al. fix the faint- 
end slope at q = —1.25 for the ellipticals, and at Q = —1.0 
for the spirals. They find that there is a slow brightening of 
the elliptical population from z = to z — 0.54 by 
AM ~ — 0.3. The spiral population does, however, not 
show any clear evolution with redshift. A change of cosmol- 
ogy affects these results somewhat. Adopting SIa/ = 0.3 
and JIa ~ 0.7, instead of ^Im = 1.0 as in Small et al., 
increases the brightening of the ellipticals to AM ~ — 0.7. 
There will in this case also be a marginal brightening for 
the spiral galaxies. A somewh at lar ger evolutionary trend is 
found by Kodama & Bower (2001), who estimate that the 



Kauffmann 1995; Diaferio et al. 2001; EUingson et al. 2001 



bright spiral population have faded by about one magnitude 
since z ~ 0.4. The fading of the bright population is there- 
fore smaller than of the faint population, which is in turn 
consistent with the disappearance of the plateau. 



5.2 The radial population gradient 

As we have just discussed, galaxies are in a hierarchical sce- 
nario accreted from the field on to the clusters (e.g., Diaferio 
et al. 2001). When field galaxies, which are predominately 
spirals with ongoing star formation, are accreted, star for- 
mation is halted after a possible burst of star formation. The 
quenching can either be abrupt, or there could be a more 
gradual decline in star formation (e.g. Barger et al. 1996, 
Poggianti et al. 1999, Balogh et al. 1999, Balogh, Navarro 
& Morris 2000, Diaferio et al. 2001). The decrease in star 
formation should set up a radial gradient in the clusters, 
where blue galaxies are expected to be more numerous in 
the outer parts, as is confirmed by observations and simula- 
tions (Diaferio et al. 2001). Also, the higher accretion rate 
at earlier epochs, as well as higher star formation rate in the 
general field, implies that the fraction of blue galaxies should 
be higher in high-z clusters. There is therefore a connection 
between the shape of the LF and the existence of a radial 
colour gradient. 

In the right panel of Fig. ^ we showed the radial distri- 
bution of late-type galaxies, divided into a bright subsample, 
with R < 23, and a faint subsample with 23 < 7i < 25. 
From the figure it is obvious that bright, blue galaxies are 
rare across the whole cluster. The flat, or at most slightly 



© 2001 RAS, MNRAS 000, 



16 Dahlen, Fransson & Ndslund 



increasing, projected surface density is consistent with a 
decline in the spatial distribution towards the core, even 
though the small statistics do not make this significant. 

The faint blue galaxies show a more dramatic trend. 
At radii ^ 0.4 Mpc this is the most numerous population 
in the cluster. Here, the surface density is consistent with 
a volume density that is either uniform, or decreasing with 
radius. The drop in surface density at the innermost point 
is consistent with an absense of faint blue galaxies at the 
core of the cluster. An explanantion may be provided by 
'galaxy harassment', which should be increasingly efficient 
in the denser core region. 

The increasing infall rate at high redshifts, predicted 
by hierachical clustering, should result also in a change with 
redshift of the radial gradient of the frac tion o f early-type 
galaxies in the cluster. Ellingson et al. (2001) show that 
galaxies with an old population dominate to larger radii in 
low redshift (0.18 < z < 0.30) clusters, giving a flatter ra- 
dial gradient, compared to high redshift (0.30 < z < 0.55) 
clusters. This indicates that high-z clusters have a higher 
population of recently accreted, star forming galaxies at 
radii closer to the core, causing a steeper radial gradient. 
As the number of recently accreted star forming galaxies 
declines with time, the gradient becomes less steep. 

With similar magnitude limits and the same cosmology 
as in Ellingson et al., we find that the radial gradient of 
C11601-f42 is at least as steep as the slope they find from 
their high redshift sample, and significantly steeper than 
that of the low redshift clusters. 



5.3 The radial and magnitude dependence of the 
Butcher-Oemler effect 

To check whether the blue fraction, /s, decreases towards 
the core as a result of a decreasing star formation rate in 
the accreted field galaxies as they fall deeper into the clus- 
ter potential, we plot in Fig. h]l fs as function of radius 
(here given in terms of R30 = 0.6 Mpc). From this figure it 
is clear that there is a trend for an increasing value of /s 
at larger radii. Th is is in agreement with both B084 and 
Kodama & Bower (2001), and consistent with the infall sce- 
nario. As is discussed by Kodama & Bower, the change of 
/b with radius further illustrates the necessity of a correct 
determination of i?3o for a meaningful comparison of the 
blue fraction between clusters. 

It is also interesting to calculate the blue fraction as 
a functi on of the magnitude cut, as is done in Kodama & 
Bower (2001) for their seven clusters at 0.23 < z < 0.43. 
When varying the magnitude cut by ± 1 mag from the fidu- 
cial value Mv = — 20 they find that /s remains almost 



constant. Only for two of the most distant clusters is there 
a trend for /s to increase at fainter magnitudes, which Ko- 
dama & Bower, however, attribute to increasing photometric 
errors with magnitude. Instead, they argue that /s should 
depend only weakly with magnitude. 

In Fig. ^we plot fs as a function of the magnitude cut. 
We find a clear trend for an increasing value of /s as fainter 
galaxies are included, similar to the trend for two of the 
highest redshift clusters in Kodama & Bower, but contra- 
dicting the general claim by Kodama & Bower that the blue 
fraction should be constant with magnitude. As illustrated 
by Fig. O, the increasing value of /s at faint magnitudes in 
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Figure 11. The blue fraction, /s, as a function of radius. R30 is 
the fiducial value used by B084 



C1160H-42 is a direct consequence of the different shapes of 
the red and blue LFs. The steepness of the blue LF causes 
the blue fraction to increase at fainter magnitude cuts, and 
decrease at bright. 

Because we are more than one magnitude brighter that 
our limit, the increase in photometric errors is <^ 0.03 mag 
over the magnitude range where we calculate fs ■ It is there- 
fore unlikely that the strong trend of /s with magnitude in 
C11601-f42 can be explained by increasing photometric er- 
rors at fainter magnitudes, as proposed by Kodama & Bower 
for the clusters showing this trend in their investigation. 

Note also that systematic effects when selecting the BO- 
galaxies may lead to an underestimate of the blue fraction 
in the faintest bin by A/s ~ 0.05, as discussed in Section 
4.4. This means that the trend shown in Figure |l^ could be 
even stronger. 

The dependence on the magnitude cut makes the cal- 
culated blue fraction dependent on the choice of cosmology. 
E.g., a constant cut in absolute magnitude, in the B084 case 
Mv = —20, corresponds to a fainter apparent magnitude in 
a A-dominated cosmology than in a cosmology without cos- 
mological constant. This should lead to a relatively higher 
blue fraction at increasing redshifts for the A cosmology. Us- 
ing Q.M = 0.3 and Q,a = 0.7, instead of the cosmology used 
by B084 {yiM ~ 0.2, SIa = 0) changes the cut in apparent 
magnitude by Am = 0.24 a,t z — 0.54. This increases the 
blue fraction from /s = 0.15 to /s — 0.18. The difference 
between the extrapolated fs sX z = 0.54 from B084, and 
the /s found for C11601-I-42 should therefore decrease in a 
A-cosmology. 

We finally notice that the blue fraction for the field, 
derived from the EIS data, is /s ~ 0.40 ± 0.07 at the 
cluster redshift (0.5 li z ^ 0.6), and /s ~ 0.39 ± 0.08 
at lower redshift [0.1 <^ z <^ 0.4). This is close to the blue 
fraction /a ~ 0.41 ± 0.10 calculated for the field by 
B084. The similarity between the blue fraction in the low 
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Figure 12. The fraction of blue galaxies, fs, for C11601+42 as a 
function of limiting magnitude. The fiducial value used by B084 
is Mv=-20. 



and intermediate redshift bins does not reflect the strong 
increase with redshift in the s tar formation in field galaxies 
found by Diaferio et al. (2001 ). This is somewhat surprising, 



even though the relatively large errors in /s do not rule out 
such trend. The fraction of late-type galaxies for the field is 
~ 0.7. This fraction is consistent with the rough distribution 
of field galaxy types used by Small et al. (1997). 



5.4 Metallicity versus star formation 

Colours of galaxies do not only refiect star formation, but are 
also dependent on metallicity, reflected in the CM relation 
(e.g. Faber 1977; Arimoto & Yoshii 1987). In principle, a low 
metallicity could cause faint old, but blue, galaxies in our 
sample to be classifled as star forming late-type galaxies. To 
estimate the importance of this effect for the blue population 
at Mb<^ — 20, we calculate the CM relation for the bright 
population with Mb < ^20, and extrapolate this to fainter 
magnitudes. 

The slope of the B — V vs. V CM relation we find for 
the early-type galaxies in C1160H~42 is A(_B - V) /AMv = 
-0.016, for Mv < - 20.7 (Mb ^ - 20.0). Over the full 
interval —22 < Mb < — 18 the mean bluening is then 
A{B — V)^ 0.07. This is small compared to the mean off- 
set in colours between early-types and late-types, which is 
A{B - 1/) ~ 0.4 at Mb ^ - 18. 

In addition, if part of the bluest old galaxies at faint 
magnitudes (—20 < AIb < — 18) are identified as late, star 
forming galaxies due to a low metallicity, the mean colour 
of the remaining early-type galaxies in this magnitude range 
will be redder than expected from an extrapolation of the 
CM relation. We do not see this effect, and therefore con- 
clude that the bluening of early-type galaxies at faint mag- 
nitudes due to metallicity should not contribute more than 



marginally to the faint blue population we detect. The blue 
colours should therefore reflect the star formation. 

Further s uppor t for this conclusion comes from Fer- 
reras & Silk (200C) who determine rest frame NUV (~ 
2000A) to V colours of early-type galaxies in C10939-I-4713 
at z = 0.42. Comparing the observed colours to synthetic 
models with different star formation histories, metallicities 
and dust models, they find that the bluening at faint magni- 
tudes is consistent with at least a fraction of the stars in the 
galaxies being young. Low metallicities alone are insufficient 
to explain the observed colours. (For further discussions on 
metallicity effects in cluster galaxies see e.g., Couch et al. 
1998; Dressier et al. 1999; Poggianti et al. 1999.) 



6 CONCLUSIONS 

Using photometry in five broad-band filters we deter- 
mine photometric redshifts for the galaxies in the field of 
C11601-I-42 at z ~ 0.54. Cluster membership can therefore 
be assigned with a higher reliability for the individual galax- 
ies, which, compared to statistical methods, significantly re- 
duces the amount of background subtraction needed. As a 
by-product we also obtain information about the spectral 
properties of the galaxies, making assumptions about e.g. 
K-corrections unnecessary. Using this method, we calculate 
the luminosity function, the radial distribution of galaxies 
within the cluster, as well as the distribution of early-type 
and late-type galaxies. With this knowledge we can directly 
understand e.g., the dependence of the Butcher-Oemler ef- 
fect on such parameters as radius and magnitude cut. 

We find a steep faint-end slope of the rest frame B- 
band cluster luminosity function, with a ~ — 1.4. This 
is consistent with the slope at faint magnitudes for both lo- 
cal and 0.1 <, z <y 0.3 clusters. The steepness of the LF is 
found to be correlated to the radial distance from the clus- 
ter centre, with a more shallow slope near the core, com- 
pared to the outer parts of the cluster. This is closely re- 
lated to the difference found in the radial distribution be- 
tween bright and faint galaxies, where bright galaxies are 
more centrally concentrated. Most surveys reporting flatter 
slopes, —1.2f^a^ — 1.0, only cover the core region of 
the cluster, or do not probe the faint population beyond the 
upturn in the LF. 

When we divide the cluster galaxies spectroscopically 
into early-type and late-type galaxies, we find that these 
populations have very different LFs. Early-type galaxies, in- 
cluding ellipticals and red spirals, are best represented by a 
Gaussian. The late-type galaxies, including blue spirals, ir- 
regulars and dwarfs, are best fitted by a Schechter function 
with a steep faint-end slope, a ~ — 1.7. 

We find that the 'plateau' region of the LF at interme- 
diate magnitudes, found in most local and low-z clusters, is 
absent in C1160H-42, and we suggest that this can be due 
to a general brightening of the blue population of galaxies, 
caused by a higher star formation rate. If the progenitors 
of todays dwarf galaxies in local clusters have faded by ~ 
2-3 magnitudes since z ~ 0.5, this class of galaxies could 
contribute to the steep slope seen in C11601-I-42. 

We find significant differences in the radial distributions 
when we divide our sample into bright and faint galaxies. 
Bright galaxies, which are predominantly early- types, are 
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highly concentrated to the centre. This shows that the high 
density of early-type galaxies in the core of local clusters is 
established already at 2 ~ 0.5. The faint galaxies have a 
more shallow radial decline at increasing radii. At the core 
the surface density flattens. 

Bright, blue, late-type galaxies are rare in the cluster. 
This is consistent with the hierarchical scenario where in- 
falling spiral galaxies have their star formation quenched as 
they approach the centre of the cluster. Faint blue galax- 
ies are more immerous, and are present across the whole 
cluster. There is, however, a decrease in the surface density 
at the core, which is consistent with a region close to the 
core depleted of faint, blue star forming galaxies. A possi- 
ble mechanism could be the destruction of small galaxies 
through 'galaxy harassment'. 

The fraction of blue galaxies in C11601-I-42 is 
/s ~ 0.15, which is higher than the local value, and qual- 
itatively in agreement with the BO effect. The fraction is, 
however, lower than predicted from the average /s(z) rela- 
tion found by B084. The blue fraction is found to depend 
on the distance to the core centre, in accordance with other 
studies. The increasing value of /s with increasing radius is 
consistent with the hierarchical scenario, where clusters are 
formed by accretion of blue field galaxies, /s is also found to 
increase with the limiting magnitude. This is a direct conse- 
quence of the different LF's for the early-type and late-type 
galaxies. The radial dependence and the shape of the LF, 
especially its faint population, are therefore connected. 

Finally, we note that although this study goes deeper 
than most previous investigations at this redshift, the con- 
clusion s are only based on one cluster, and a large sample 



including clusters of varying richness and concentration is 



needed to access the robustness of our conclusions. We hope 
to address this in future studies. 
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